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, We compute opacities for the electronic molecular band systems A ^11+ - 

^ ' X of CrH and CrD, and A"^!!- X of MgH and MgD. The opacities 

are computed by making use of existing spectroscopic constants for MgH and 
CrH. These constants are adjusted for the different reduced masses of MgD 
, and CrD. Frank-Condon factors are used to provide intensities for the indi- 

. vidual vibronic bands. These results are used in the computation of synthetic 

spectra between Toff — 1800 and 1200 K with an emphasis on the realisa- 
tion of "deuterium test", first proposed by Bejar et al. (1999) to distinguish 
brown dwarfs from planetary mass objects. We discuss the possible use of CrD 
and MgD electronic bands for the "deuterium test" . We find CrD to be the 
more promising of the two deuterides, potentially, the most useful bands of 
CrH/CrD are the Av = +1 and Aw = -1 at 0.795 and 0.968 /im . 
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- planets: evolution 



* E-mail:yp@mao. kiev.ua 



2 Ya. V. Pavlenko et al. 
1 INTRODUCTION 

The "deuterium test" was suggested as a method of identifying planetary mass objects 
among cool objects (Bejar et al. 1999, Chabrier et al. 2000). In practise, it was proposed 
to search for absorption lines of molecules containing deuterium (HDD, CrD, FeD, etc.). 
Deuterium is burnt in stellar interiors via the fusion reaction ^D(p,7)^He at temperatures 
(T > 8 X 10^ K). The interiors of substellar objects with M < 13 Mj , where Mj is a Jovian 
mass (0.001 Mq ), do not reach temperatures high enough for deuterium to ignite (Saumon 
et al. 1996). As a result the deuterium abundance in the atmospheres of these objects is 
unchanged from the formation of these objects. This gives rise to the definition of a brown 
dwarf as an object which has insufficient mass to fuse to ^He, but has sufficient mass 
to fuse to ^He. By comparison, a planetary mass object has insufficient mass to ignite 
fusion of any sort (Saumon et al. 1996). 

In higher-mass objects such as stars, deuterium burning is completed comparatively 
quickly (t = 1 - 3 million years) during the evolution prior to the star's period on the main 
sequence (D'Antona & MaziteUi 1998). The deuterium depletion rate depends on the mass 
of the star or brown dwarf and thus the "deuterium test" can be used in a number of ways. 

• Discern planets from brown dwarfs in a population of low-mass objects. 

• Determine the evolutionary status of young objects in open clusters with ages of several 
million years. 

• Study the evolution of the abundance of deuterium in the atmospheres of low-mass 
substellar objects, a phenomenon which is poorly understood. The rate of depletion of deu- 
terium depends upon rotation, magnetic field strength, and other parameters which affect 
the efficiency of convection in low mass objects. 

Such investigations can usefully be combined with the "lithium test" proposed by Rebolo 
et al. (1992), Magazzii et al. (1993). For stars (M > 75 Mj ), the burning of lithium, Li 
(p,Q;)'^He, becomes efficient at evolutionary stages preceding the main sequence at interior 
temperatures of T > 2.5 x 10^ K (D'Antona & Mazitelh 1998). The "lithium test" has been 
successfully applied to identify brown dwarfs in a population of ultracool dwarf stars (Rebolo 
et al. 1996, Basri 2000). The lithium test is relatively easily applied to M-dwarfs. The reason 
for this is that the resonance lines of neutral atomic lithium lie in the optical region of 
the spectrum, at 0.6708 fim. . One blemish with the "lithium test" is the severe blending 
of lithium lines with the background of TiO lines, nonetheless Li lines break through the 
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molecular background. In the spectra of cooler M and L dwarfs Ti atoms are depleted onto 
dust particles (Tsuji et al. 1996, Jones & Tsuji 1997, Pavlenko 1998) so TiO absorption is 
weakened. In the L dwarf regime the appearance of lithium lines contends with dust opacity 
and the wings of K I and Na I resonance lines (see Pavlenko et al. 2000 and references 
therein) . 

The realisation of the "deuterium test" is considerably more difficult. Conventionally, 
the deuterium abundance of hot objects is determined from analysis of multicomponent 
features on the background of the 0.1215 ixm H I line seen in emission. Unfortunately, 
this method cannot be used in the case of ultracool dwarfs which are covered by a thick 
envelope of neutral hydrogen. Observations also yield H I emission hues {Ha 0.6563 jjim ) 
formed in the outermost layers of hot chromospheres or accretion disks around young stars. 
However, such a plasma is variable and may be polluted by interstellar material. 

Taking into account these circumstances, it is logical to analyse the spectral lines of 
deuterated molecules formed in the photospheric layers of ultracool dwarfs. The first inves- 
tigation and analysis of the combined spectra of H2O /HDO were carried out by Chabrier et 
al. (2000) and Pavlenko (2002). Due to the change in mass and the breakdown of molecular 
symmetry the vibration-rotation bands of HDO in the mid-infrared spectrum shift with re- 
spect to the H2O bands. There are a few spectral regions which can be used for detection of 
HDO lines in the IR spectra of ultracool dwarfs: 3.5 - 4, ~5, 6-7 ^m. (Pavlenko 2002). The 
main problem is that despite the shift in wavelength such HDO lines will be on a background 
of far stronger H2O lines. 

A possible alternative to HDO are the diatomic hydrides. Strong molecular bands of 
diatomic metal hydrides such as MgH and CrH can be observed in the optical spectrum of 
ultracool dwarf stars. The MgH band system A'^H - X ^E"*" can be observed at 0.47-0.6 
/xm , and the CrH band system A - X ®S+ show absorption features at 0.6-1.5 //m . 

The molecules MgH and CrH have been known in astrophysics for a long time. MgH 
has been more extensively studied than CrH, because it can be observed in the spectra of 
G to M type stars. The dissociation energy of MgH is very low (1.285 eV) so lines of this 
molecule are very sensitive to the temperature and gravity variations in stellar atmospheres. 
MgH lines were used to determine temperatures in the atmospheres of cool giants (Wyller 
1961) and the Sun (Sinha et al. 1979, Sinha & Joshi 1982), and for the determination of the 
surface gravity of stars (Bell & Gustaffson 1981, Bell at al. 1985, Berdyugina & Savanov 
1992, Bonnel & Bell 1993). 
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The pure rotational spectrum of MgH and MgD radicals (X ^S"^ ) in their ground state 
^;=0 and v=l vibrational modes has been studied by Ziurys et al. (1993). The first MgH 
linehst was computed by Kurucz (1993). Recently, more extensive studies of MgH transitions 
were performed by Week et al. (2003, 2003a, 2003b) and Skory et al. (2003). 

Although CrH has been known since Gay don & Pearse (1937), its electronic spectrum 
remained relatively unstudied for many years. Engvold ct al. (1980) identified lines of CrH 
in a spectrum of sunspots as formed by X - X transitions. They used the results of 
studies of multiplicity of E-terms of CrH by Klehman & Uhler (1959) and O'Connor (1967). 
Later Ram, Jarman & Bernath (1993) performed a rotational analysis of 0-0 band of the 
A - X electronic transition and obtained improved rotational constants for the v' 
=0 vibrational state. Combining these results with those of Bauschlicher et al. (2001) and 
Lipus et al. (1991) for the vibrationally excited transitions. Burrows et al.(2002) computed 
an extended linelist for CrH. 

Recently Shin et al. (2005) have measured radiative hfetimes of the v — 0,1 levels of 
A ^E+ state of CrH. These measured lifetimes are about 16% - 45 % longer that those ob- 
tained by Burrows at al. (2002). These results provide evidence that the oscillator strengths 
of Burrows et al. (2002) should be corrected by a factor of 0.8 for at least the transitions to 
the v' = 0. 

The submillimeter spectra of CrH and CrD formed by pure rotational transitions in the 
ground electronic state, have been observed in the laboratory by Halfen & Ziurys (2004). 
Electronic bands of MgD and CrD are likely to be located in the same spectral regions 

as the corresponding bands of MgH and CrH. In this paper we model the bands of these 
molecules to analyse the possibility of their use for the determination of the D/H ratio in 
the atmospheres of ultracool dwarfs. 

In section 2 we present a description of the procedures used to compute the molecular 
bands of CrH, CrD, MgH and MgD. In section 3 we present the vibrational- rotational 

constants of MgH, MgD, CrH and CrD. In section 4 we present the results of the computation 
of molecular bands. In section 5 we discuss the possibility of using the electronic bands of 
diatomic molecules for a realisation of the "deuterium test" . 
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2 PROCEDURE OF COMPUTATIONS 



2.1 Profile of electron bands of diatomic molecule 

In this paper we compute the profile of the electronic band averaged over its rotational 
structure (Naersisian et al. 1989): 



where as symbols " label lower levels of transitions, A — S7r^{3hc)^^, k is Boltzmann's 
constant and T is the temperature, a; is a frequency and are Pranck-Condon factors for 
the corresponding transition, the strength of transition — 5'e(0, 0) x a;,;'„"/<^o,o (Schadee 
1968), Q{T) is the partition function of the molecule. 

For a given frequency the sum of the contribution from the P, Q and R branches is: 



where Fe = exp{-hcE'^ /{kT))/5Ee, 

C,j = + 5E,)(2j" + 1) exp(-/ic(4' + E'^)), 

X^l/Su. 

This method of the computation of the profiles of molecular bands was successfully 
applied to the modelling of the bands of TiO in the spectra of ultracool dwarfs (Pavlenko 
1997a, 1998a, Pavlenko 2000), and CN and C2 bands in the spectra of evolved stars (Pavlenko 
et al. 2000a,c). It is worth noting that direct comparison of our JOLA synthetic spectra with 
more sophisticated "line by line" computations for linelists of a few molecules show good 
agreement (Pavlenko 1998a) both in the profiles of bands and positions of the band heads. 

The JOLA approach has been used for the modeUing of the CrH bands Av =0 band sys- 
tem located at 0.8640 fim in the spectra of ultracool dwarfs (Pavlenko 1999, Pavlenko et al. 
2000). These bands are not considered for the deuterium test here, because the wavelengths 
of corresponding lines of CrH and CrD are nearly identical. 




(1) 



ASeFe 

n\Q{T) 




(2) 
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2.2 Computations of molecular constants of isotopic molecules 

For a given electronic state, the energy levels of a diatomic molecule can, in general, be fitted 
by: 

Ee,v,N = T^ + U;e(v + l/2)-U;eXe{v + l/2f 

+BeN{N + 1) - aeN{N + l){v + 1/2) 
-DeiN{N+l)f 

+(3,{N{N + l))\v + l/2), (3) 

where Te, Ue, oJeXe, -Be, <^e, and Pe are the various spectroscopic constants, v and N are 
vibrational excitation and rotational quantum numbers, respectively. To a good approxima- 
tion the electronic structure of the various isotopic analogues of a given molecule are the 
same. However, as the nuclei have different masses, the reduced mass and moment of inertia 
of the molecule are different, resulting in different rotation-vibration frequencies. For the 
case of an anharmonic oscillator and non-rigid rotator model of the molecule, the rotational- 
vibrational constants depend upon the ratio of reduced mass fiab = iTia x mb/{ma + nib) of 
two isotopic molecules (see Wang & Xia 1996, Craybeal 1988). 

The ratios of spectroscopic constants of two isotopic species, labelled as "o" and "i", 
of reduced masses fXo and /Xj, depend on the ratio p — {jio/ l^i)^^^'^^ ■ So that the following 
approximations can be used: 

{uJeXe)i/{uJeXe)o = 



(ae)i/(«e)o = P^ 
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Table 1. Reduced masses and ratios of reduced masses for molecules of interest. 



Molecule 



P 



MgH 



0.967481 



0.721116 



MgD 



1.86051 



CrH 



0.988438 



0.713873 



CrD 



1.93958 



m^/{Pe)o = P' 



The reduced masses of the molecules of interest here are given in Table [T] using the 
dominant species ^^Mg and ^^Cr. The spectroscopic constants of MgH and CrH have been 
taken from Balfour & Cartwright (1976), Lemoine et al (1988) and Bauschlicher et al. (2001). 
The computed spectroscopic constants for MgD and CrD are given in Table [21 We adopted 
oscillators strengths fg = 0.059 (Kuznezova et al. 1980) and 0.001 (Pavlenko 1999) for the 
MgH {A^U- X ) and the CrH {A - X ) band system, respectively. 

We calculated the Franck-Condon factors q^/^/r for the MgH and CrH band profile com- 
putations. The RADEN program of Kuzmenko et al. (1984), and the FRANKQ program 
(Tsymbal, 1984) were used, respectively. In the FRANKQ program a Morse potential is used 
to represent the true potential, and is solved for rotation-vibration motion. We compared the 
computed values with the results of Nicholls (1981) and found only small differences. 
It is worth noting that the adopted system of determines the relative strength of the 
(f ",f ') bands, but not their location in the spectrum. In this work we are interested in the 
computation of the main bands of molecules of interest, therefore the choice of any g^/^,// does 
not critically affect our main results. 

2.3 Partition function and equilibrium constants 

The internal partition function is given by a sum over all states: 



where, S is electronic spin, Ee,v,N is the energy of the state with electronic excitation e, 
vibrational excitation v and with rotational angular momentum N, Eq is the zero point 
energy. When rotational angular momentum is coupled to the electronic spin a splitting of 




(4) 
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Table 2. Spectroscopic constants of MgH, MgD, CrH and CrD in cm~^. 



Molecule, 



data source 












Be 








/5e 




MgH, 
























Balfour & Cartwright (1976) 


A 


2H 


19216.8 


1599.50 


32.536 


6.1913 


0.1931 


3.60x10- 


-4 


6.1x10" 


6 


Lemoine et al. (1988) 


X 


^2+ 


0.00 


1495.25 


31.637 


5.6443 


0.1845 


3.53x10" 




2.57x10 


— 6 


MgD, 
























This work 


A 


2H 


19216.8 


1153.43 


16.918 


3.2195 


0.07241 


9.73x10" 


-5 


1.2x10" 


6 


Lemoine et al. (1988) 


X 


2S + 


0.0 


1078.14 


16.147 


2.9668 


0.068294 


9.61x10" 


-5 


2.79x10 


-7 


CrH, 
























Bauschlicher et al. (2001) 


A 




11552.7 


1524.80 


22.280 


5.3427 


0.141349 


3.01x10" 


-4 


9.18x10- 


-5 




X 


6S+ 


0.0 


1656.05 


30.491 


6.2222 


0.180978 


3.52x10" 


-4 


5.21x10" 


-6 


CrD, 
























this work 


A 


6S+ 


11552.7 


1088.32 


11.35 


2.7227 


0.051395 


7.82x10" 


-5 


1.70x10 


-5 




X 


6E+ 


0.0 


1182.00 


15.533 


3.1709 


0.065804 


9.14x10" 


-5 


9.66x10 


-7 


Ram k. Bernath (1995) 


A 


6E+ 


11559.7 


1066.42 




2.7627 


0.049687 


7.03x10" 


-5 








X 


6E+ 


11559.7 


1183.19 


15.60 


3.1754 


0.065908 


9.09x10" 


-5 







the spin degeneracy occurs. For the case of the CrH X electronic states, S=5/2, this 
results in 6 separate states with angular momentum J = N±5/2, N±3/2, N±l/2. Similarly 
for MgH, X J = N±l/2. The internal partition function is then expressed as: 

Q.., = E (2^ + 1) exp (z(E,-j^JlM\ (5) 

e,v,J,N \ / 

Energy levels are computed using the rotation, vibration constants listed in Table [2l It is 
worth noting a good agreement of our data with the constants obtained by Ram & Bernath 
(1995) from the high resolution spectroscopy of the A ^S"*" - X system of CrD molecule. 
The partition functions of MgH, MgD, CrD are computed via direct summation of these 
levels, see equation |H For CrH, we have obtained accurate energy levels from the linelist 
of Burrows et al. (2002). This list of energy levels covers the X ^E"*" and A ^E"*" electronic 
states, extends from J = 0.5 to 39.5 and includes the vibrational states v =0,1,2,3 in the 
ground electronic state and v =0,1,2 for the excited electronic state. The CrH partition 
function is calculated by the direct summation of the Burrows et al. (2002) energy levels. In 
addition, to cover the states omitted from the Burrows et al. (2002) linelist, we supplement 
energy levels computed with the spectroscopic constants. The partition functions are listed 
in Table [31 

Figures [U and [2] show the partition functions and the respective fits of Sauval & Tatum 
(1984). For CrH there is excellent agreement with the partition function of Sauval & Tatum 
(1984), but for MgH our partition function is slightly larger than that of Sauval & Tatum 
(1984) at low temperature and smaller at temperatures above 3200 K. Sauval & Tatum 
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Table 3. Internal partition functions, Q, as a function of temperature (powers of ten are given in parenthesis). 



T/K Q(MgH) Q(MgD) Q(CrH) Q(CrD) 



50 


1 


322 


1) 


2.440 


1) 


3.608 


1) 


6.855 


1) 


100 


2 


575 


1) 


4.815 


1) 


7.011 


1) 


1.351 


[2) 


150 


3 


832 


1) 


7.194 


1) 


1.042 


2) 


2.018 


[2) 


200 


5 


091 


1) 


9.582 


1) 


1.384 


2) 


2.687 


'2) 


250 


6 


353 


1) 


1.199 


2) 


1.726 


2) 


3.359 


2) 


300 


7 


623 


1) 


1.445 


2) 


2.069 


2) 


4.041 


2) 


350 


8 


905 


1) 


1.698 


2) 


2.415 


2) 


4.738 


2) 


400 


1 


021 


2) 


1.961 


2) 


2.764 


2) 


5.455 


2) 


450 


1 


154 


2) 


2.236 


2) 


3.118 


2) 


6.199 


2) 


500 


1 


290 


2) 


2.524 


2) 


3.479 


2) 


6.973 


2) 


550 


1 


431 


2) 


2.825 


2) 


3.849 


2) 


7.780 


[2) 


600 


1 


577 


2) 


3.142 


2) 


4.228 


2) 


8.623 


[2) 


650 


1 


727 


2) 


3.474 


2) 


4.618 


2) 


9.503 


'2) 


700 


1 


883 


2) 


3.823 


2) 


5.020 


2) 


1.042 


3) 


750 


2 


045 


2) 


4.188 


2) 


5.435 


2) 


1.138 


3) 


800 


2 


212 


2) 


4.569 


2) 


5.864 


2) 


1.238 


3) 


850 


2 


386 


2) 


4.969 


2) 


6.306 


2) 


1.343 


3) 


900 


2 


566 


2) 


5.385 


2) 


6.763 


2) 


1.451 


3) 


950 


2 


753 


2) 


5.820 


2) 


7.235 


2) 


1.564 


3) 


1000 


2 


947 


2) 


6.272 


2) 


7.722 


2) 


1.682 


.3) 


1100 


3 


355 


2) 


7.233 


2) 


8.744 


2) 


1.930 


;3) 


1200 


3 


791 


2) 


8.267 


2) 


9.831 


2) 


2.197 


'3) 


1300 


4 


256 


2) 


9.378 


2) 


1.098 


3) 


2.482 


3) 


1400 


4 


750 


2) 


1.057 


3) 


1.221 


3) 


2.787 


3) 


1500 


5 


274 


2) 


1.183 


3) 


1.350 


3) 


3.111 


3) 


1600 


5 


827 


2) 


1.318 


3) 


1.487 


3) 


3.455 


3) 


1700 


6 


410 


2) 


1.460 


3) 


1.630 


3) 


3.818 


3) 


1800 


7 


023 


2) 


1.610 


3) 


1.782 


3) 


4.203 


3) 


1900 


7 


664 


2) 


1.767 


3) 


1.941 


3) 


4.608 


.3) 


2000 


8 


333 


2) 


1.932 


3) 


2.108 


3) 


5.035 


;3) 


2100 


9 


029 


2) 


2.104 


3) 


2.283 


3) 


5.483 


'3) 


2200 


9 


751 


2) 


2.283 


3) 


2.467 


3) 


5.953 


3) 


2300 


1 


050 


3) 


2.469 


3) 


2.659 


3) 


6.446 


3) 


2400 


1 


127 


3) 


2.661 


3) 


2.859 


3) 


6.962 


3) 


2500 


1 


206 


3) 


2.859 


3) 


3.069 


3) 


7.501 


3) 


2600 


1 


287 


3) 


3.063 


3) 


3.288 


3) 


8.063 


3) 


2700 


1 


370 


3) 


3.272 


3) 


3.517 


3) 


8.650 


3) 


2800 


1 


455 


3) 


3.486 


3) 


3.755 


3) 


9.261 


.3) 


2900 


1 


542 


3) 


3.705 


3) 


4.003 


3) 


9.896 


;3) 


3000 


1 


629 


3) 


3.928 


3) 


4.261 


3) 


1.056 


'4) 


3100 


1 


719 


3) 


4.154 


3) 


4.530 


3) 


1.124 


'4) 


3200 


1 


809 


3) 


4.385 


3) 


4.810 


3) 


1.195 


4) 


3300 


1 


901 


3) 


4.618 


3) 


5.101 


3) 


1.269 


4) 


3400 


1 


993 


3) 


4.854 


3) 


5.403 


3) 


1.345 


4) 


3500 


2 


086 


3) 


5.093 


3) 


5.716 


3) 


1.424 


4) 


3600 


2 


180 


3) 


5.334 


3) 


6.042 


3) 


1.506 


4) 


3700 


2 


275 


3) 


5.577 


3) 


6.379 


3) 


1.590 


4) 


3800 


2 


370 


3) 


5.821 


3) 


6.728 


3) 


1.676 


;4) 


3900 


2 


465 


3) 


6.067 


3) 


7.089 


3) 


1.765 


;4) 


4000 


2 


561 


3) 


6.314 


3) 


7.463 


3) 


1.857 


;4) 


4100 


2 


657 


3) 


6.562 


3) 


7.849 


3) 


1.951 


'4) 


4200 


2 


753 


3) 


6.811 


3) 


8.248 


3) 


2.048 


'4) 


4300 


2 


849 


3) 


7.060 


3) 


8.659 


3) 


2.148 


'4) 


4400 


2 


945 


3) 


7.310 


3) 


9.084 


3) 


2.250 


'4) 


4500 


3 


041 


3) 


7.559 


3) 


9.521 


3) 


2.354 


4) 


4600 


3 


137 


3) 


7.809 


3) 


9.971 


3) 


2.461 


4) 


4700 


3 


233 


3) 


8.059 


3) 


1.043 


4) 


2.571 


4) 


4800 


3 


329 


3) 


8.308 


3) 


1.091 


4) 


2.683 


'A 


4900 


3 


424 


3) 


8.557 


3) 


1.140 


4) 


2.797 


'A 


5000 


3 


519 


3) 


8.805 


3) 


1.190 


4) 


2.913 


;4) 
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Figure 1. Partition functions of MgH and MgD computed in tiiis work, and the fit of Sauval & Tatum (1984). 



(1984) make use of a high temperature approximation to the rotational partition function. 
This uses the rigid rotor approximation and is determined by rewriting the partition function 
sum as an integral from zero to infinity. In contrast our partition function sum is truncated 
at the dissociation energy Dq, which for MgH is only 1.285 eV. Thus at high temperature 
Sauval & Tatum (1984) predict a higher partition function than would be expected by direct 
summation of bound states. 

We have calculated equilibrium constants for the reactions: MgH+D^MgD+H and 
CrH+D^CrD+H. The dimensionless equilibrium constants are given by: 
n{CrH)n{D) 



K 



n{CrD)n{H) 
Q{CrH)Q{D) 



exp 



(6) 



Q{CrD)Q{H) 

where n are the number densities, and AE is the energy difference between zero point 
energies of CrH and CrD; Q is the total partition functions which can be factorised into 
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Figure 2. Partition functions of CrH and CrD computed in tliis work, and tiie fit of of Sauval and Tatum (1984). 

internal and translational components, Q = QintQtran, where Qtran = (27rMjA;T//i^)^/^. 
Qj is the degeneracy of energy level j. By convention nuclear spin degeneracy is omitted 
but electronic spin degeneracy is included. Thus for a proton Qint = 1 and for an electron 
Qint = 2. The equilibrium constants are listed in Table HI The equilibrium constants indicate 
that at low temperatures the formation of the deuterides are strongly favoured over the 
hydrides, but at high temperatures the hydrides become weakly favoured. 

It is worth noting that in cool dwarfs hydrogen and deuterium are primarily found 
in their molecular forms, H2, HD and D2. In a full chemical equilibrium calculation the 
partition functions and dissociation energies of H2, HD and D2 will have an affect upon the 
abundances of the molecules of interest. So a more appropriate reaction may be MgH + HD 
^ MgD + H2. However, to calculate equilibrium constant for this reaction would require 
the partition functions of Q(HD) and Q(H2) and their dissociation energies. Expressing the 
equilibrium constants relative to atomic D and H avoids any uncertainties arising from the 
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partition functions and dissociation energies of H2, HD and D2. A calculation of the partition 
functions and dissociation energies of H2, HD and D2 is beyond the scope of this project. 

In calculating the equilibrium constants we assumed that the electronic partition func- 
tions of H and D are identical. This is a good approximation as the main differences in the 
energy levels are proportional to the change in reduced mass and are thus small. 



2.4 lonisation-dissociation equilibrium 

The equations of ionisation-dissociation equilibrium were solved for media consisting of 
atoms, ions and molecules. We took into account ~ 100 components (Pavlenko 1998, 2000). 
The constants for equations of chemical balance were taken from Tsuji (1973) and Gurvitz 
et al. (1989). To find the densities of species we apply the technique applied by Kurucz 
(1975) in the ATLAS5 - ATLAS12 programs. In the framework of the local thermodynamic 
equilibrium (LTE) approach, the densities of a molecule consisting of atoms labelled 1,2, ...n 
in the k+ ionisation stage can be described by the Gouldberg-Waage equation: 

77^], (7) 



U[n, UIQ^ kT 
where is the number density of the i^^ component, Qi are the total partition functions for 
species i, and -Ei,2,...« is the dissociation energy (Dq) of molecule which consists of / atoms 
and has p"*" charge. 

In the ATLAS program a computationally more convenient format is used: 



' - exp[-E,,,,,JkT + g{T)], (8) 



A^l,2,..i 

here the function g{T) is given by: 



giT) = h~c{T + d{T-e{T + fT))) (9) 
+ ^(/-p-l)lnT, 

values of h, c, d, e (Table [5]) are determined by the fitting to the data (see Table Hj) using the 
least square minimisation procedure for the temperature range 400 - 5000 K. 

As was expected, at LTE the equilibrium constants of our hydrides and deuterides are 
very similar. Indeed, the equilibrium constants are ~ Qint,i/Mi, Mi is the mass of the i^^ 
species. Partition functions of the deuterated molecules are larger by approximately a factor 



Bands of CrD and MgD and the "deuterium test 

Table 4. Equilibrium constants for deuterium exchange reactions (powers of ten are given in parenthesis). 
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Table 5. Chemical equilibrium constants for MgH, MgD, CrH and CrD (powers of ten in parenthesis). 



Molecule Dq /eV beds f 

CrH 2.863 3.690(1) 6.678(-3) 2.541(-6) 5.115(-10) 3.929(-14) 

MgH 1.285 3.605(1) 6.834(-3) 2.551(-6) 5.049(-10) 3.834(-14) 

CrD 2.888 3.731(1) 6.919(-3) 2.607(-6) 5.192(-10) 3.967(-14) 

MgD 1.311 3.648(1) 7.111(-3) 2.647(-6) 5.214(-10) 3.942(-14) 

From Shayesteh et al. (in prep.) 

2, which corresponds to the larger molecular weight of deuterium. Fitted chemical equilib- 
rium constants computed for CrH, CrD, MgD and MgH are shown in Tabled 

2.5 Synthetic spectra 

Synthetic spectra were computed using the program WITA6 (Pavlenko 2000) assuming 
LTE, hydrostatic equilibrium and a one- dimensional model atmosphere without sources and 
sinks of energy. Hereafter we use the term "synthetic spectra" instead of "spectral energy 
distributions" to simplify the text. Theoretical synthetic spectra were computed for the 
model atmospheres of dwarfs with effective temperatures T^s = 1200 - 1800 K from the 
COND grid of Allard at al.(2001) using solar metallicity (Anders & Grevesse 1989). Unless 
otherwise mentioned all models are for log g = 5.0. 

2.6 Atomic lines opacity 

To compute the band profiles of the electronic bands in the spectra of late spectral class 
dwarfs one must account for absorption of the resonance lines of Na I (0.5891, 0.5897 /im ) 
and K I (0.7667, 0.7701 /im ), which are extremely strong. In fact, optical and near infrared 
spectra of L-and T-dwarfs are governed by them (Pavlenko 1997, 1998). Temperatures in 
the atmospheres of ultracool dwarfs are lower than 2000 K, therefore alkali atoms exist there 
mainly in the form of neutrals. Furthermore, due to the low opacity of the atmospheres of 
ultracool dwarfs their photospheres occur in high pressure layers. Strong resonance lines 
of the most abundant alkali metals, i.e sodium and potassium, are formed in very dense 
and cool plasmas. Due to the high efficiency of pressure broadening, their formally com- 
puted equivalent widths may reach a few thousand Angstroms (Pavlenko 2000). Burrows & 
Volobuev (2003) show that we cannot use the conventional theory of collisional broadening. 
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Figure 3. Molecular bands of the A ^11 — X system of MgH and MgD in the theoretical spectrum computed for a 

1800/5.0/0 COND model atmosphere. Spectra are shown for D/H = 1. The lower line on the plot shows the pure spectrum of 
MgH computed in the framework of JOLA approximation. For comparison the observed spectra of two L dwarfs 2MASS1632+19 
(Martin et al. 1999) and SDSS 0236+0048 (Leggett et al. 2001) are shown as solid lines. 



A more sophisticated theory involving quasi- stationary broadening should be used for them 
(see also Allard et al. 2003, and Zhu et al. 2006). 

In this work we use the potentials of quasi-stationary chemical interactions of K and 
Na with the most numerous species, atomic He and molecular H2 computed by GAMESS 
(Granovsky et al. 2004). Our procedure is described in more detail in Pavlenko et al. (2007). 
In calculations of K I profiles we used a combined profile: the cores of these lines were 
computed in the framework of the collisional approach and their wings {5\ > 0.004 fim ) 
were treated by quasi- stationary theory. 
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3 RESULTS 

3.1 MgH & MgD 

Computed JOLA profiles of tlie electronic band system A^Il- X of MgH and MgD are 
shown in Fig. [3l It is worth noting that a simulation based on a sound quantum mechanical 
model, including spin-orbit coupling in the A ^11 state, agreed reasonably well with the 
JOLA calculation (see Hill 2007). First, the synthetic spectra were computed for the case of 
a solar abundances (N(H)/N(total)= 0.9). Then all the hydrogen was replaced by deuterium 
(N(D)/N(Total)=0.9). As stated above, the aim of this paper is a comparative analyses of 
the possibilities of detection of the deuterated molecules MgD and CrD which are formed on 
the background of the electronic bands of the more abundant conventional hydrides MgH 
and CrH. The detection of MgD bands is very challenging for several reasons: 

• The bands of MgD lie at relatively short wavelengths. High precision observations of 
this spectral region are difficult in objects with low Tgfj . 

• The differences in the positions of the band heads of MgH and MgD are rather small. 
The MgD bands will be formed on a background of far stronger MgH bands. 

• MgH and MgD bands in the spectra of cool dwarfs are close to a strong Na I resonance 
line and the lines of other molecules/atoms. 

3.2 CrH & CrD 

The relative behaviour of CrH and CrD bands is shown in Fig. |H In general, the bands of 
the A ^S"*" - X system occupy a larger wavelength range. Due to the larger reduced 
mass, distances between heads of CrD bands of different Af are smaller than those for the 
CrH. As a result, bands of negative and positive Af are shifted bluewards and redwards 
with respect to the corresponding CrH bands. 

The transitions with Af = 1 form a band head at 0.76843 /im which coincides with the 
core of the strong K I resonance doublet. Thus in L dwarfs with strong CrH absorption the 
contribution of CrH has to be accounted in the modelling of 0.77 ^m. K I feature. 

For the Af = case the bands of CrH and CrD are co- incident. For the Ai; positive case 
the CrD band heads are displaced to longer wavelengths of the CrH band head, thus the 
strongest CrD lines will be blended with the lines from the tail of the corresponding CrH 
band. Bands with Af = 2 are blended with the K I resonance line. Bands with Af = 3 are 
located in the spectral region of the so-called "satellite rainbows" of the K I line (Burrows 
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Figure 4. Molecular bands of the A ^S"*" — X ®S+ system of CrH and CrD and a theoretical spectrum computed for a 
1800/5.0/0 COND model atmosphere. Spectra are shown for D/H =1. The lower line on the plot shows the pure spectrum of 
MgH computed in the framework of JOLA approximation. For comparison the observed spectra of two L dwarfs 2MASS1632+19 
(Martin et al. 1999) and SDSS 0236+0048 (Leggett et al. 2001) are shown as solid lines. 



& Volobuev 2003). For the Af negative case the CrD band heads are displaced to shorterer 
wavelengths of the CrH band head, these bands offer the best possibility of detection. 



4 DISCUSSION 

The detection of deuterated molecules in the spectra of ultracool dwarfs provide a challenge 
for both theoreticians and observers. Indeed, in the atmospheres of planetary mass objects 
(M < 13 Mj) we cannot expect ratio D/H > 2 x 10"^ This means the lines of deuterated 
molecules should be about 5000 times weaker than those of the hydrides. The ideal case would 
be a spectral region where the molecular bands are not blended. So a crucial requirement is 
a large difference in the wavelengths of the band heads of the hydrides and deuterides. 
CrH appears to be more useful than MgH in the search for deuterated species. The band 
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heads of CrD arc displaced significantly from the band heads of CrH. Bands of CrH are 
observed in the spectra of the latest L dwarfs (Kirkpatrick et al. 1999). The CrD bands are 
located in the "near infrared" spectra, where fluxes are much higher than in the "optical" 
spectral regions. 

In this paper we show that the most useful bands for the realisation of the deuterium 
test arc Av = +1 and Av = — 1 (A = 0.795 and 0.968 /xm , respectively). The Av = — 1 
band looks especially promising. It is located in the near infrared region with an absence of 
strong background absorption features. However, portions of this CrD band will be swamped 
by the 0-0 FeH Wing-Ford band at 0.99//m and possibly by water bands. Still, the case for 
CrD looks better than for HDO lines which are formed on a background of strong H20 lines 
(Pavlenko 2002). High quality line lists are required to test these possibilities fully. 

It is worth noting, that a potential problem lies in the possibility that Cr atoms are 
absorbed onto dust particles. The depletion of Cr will reduce the strength of both CrD and 
CrH bands. Fortunately as CrH bands are located in the same spectral region, we can "scale" 
the CrD depletion processes by adjusting the strength of CrH bands. 

For more precise studies, more accurate and detailed hnelists of CrH and CrD are re- 
quired. The calculation of such linelists would require new improved computations supported 
by new laboratory measurements. One problem that concerns us is that even once we have 
a good agreement with the model and experimental data for the MgH or CrH molecule, 
there are still perturbations about which we know very little from experiments done so far. 
Nonetheless, MgH is a non-starter for the deuterium test and as we note above, the model 
adopted in our paper is more likely to be reliable for CrH and CrD. It is worth noting 
that the use of the pure rotational-vibrational bands located in the mid and far infrared 
spectral region may offer an alternative. Indeed, the displacement between CrH and CrD 
rotation-vibration bands is even larger, than for the case of electronic bands. 

Nevertheless, we cannot be certain that we have identified the best candidate systems 
for the deuterium test. Future investigations of deuterated molecules in different spectral 
regions are important to determine which offers the best possibility for the realisation of 
this test. The ideal solution would be to detect lines of deuterated molecule(s) in different 
spectral regions. This presents a serious observational challenge which can only be met 
in combination with careful laboratory measurement and the computation of high quality 
molecular spectra. 



Bands of CrD and MgD and the "deuterium test" 19 

5 ACKNOWLEDGEMENTS 

We thank anonymous Referee for the helpful remarks. We thank Prof. Peter Bernath for a 
constructive discussion and for providing us with a copy of his paper prior to publication. 
The UK Particle Physics and Astronomy research council is thanked for visitor and PDRA 
support. The Royal Society is thanked for a visitor grant. YP is supported by a Leverhulme 
Trust grant. This work was partially supported by the Cosmo mi crophysics program of Na- 
tional Academy of Sciences and National Space Agency of Ukraine. Our research has made 
use of the SIMBAD database operated at CDS, Strasbourg, France. 



REFERENCES 

Allard, F., Hauschildt, P. H., Alexander, D. R., Tamanai, A., Schweitzer, A. 2001, ApJ, 
556, 357. 

Allard, M.F., Allard, F, Hauschildt, P.H., Kielkopf, J.F., and Machin, L. 2003, A&A., 411, 
L473. 

Anders, E., Grevesse, N., 1989, GeCoAA, 53, 197. 

DAntona, F., Mazzitelli, I. 1998. In Brown Dwarfs and Extrasolar Planets, Ed. by R. 
Rebolo, E. Marti n, and M. R. Zapatero Osorio, Astron. Soc. Pac. Conf. Ser. 134, 442. 
Balfour, W. J., Cartwright, H. M., 1976, A&AS, 26, 389 

Bauschhcher, C.W.Jr., Ram, R.S., Bernath, P.F., et al. 2001. J. Chem. Phys., 115, 1312. 
Basri, G. 2000, ARA&A, 38, 485 

Bejar, V. J. S., Zapatero Osorio, M. R., and Rebolo, R. 1999. ApJ, 521, 671. 

Bell, R. A., Gustafsson, B., 1981, BAAS, 13, 829. 

Bell, R. A., Edvardsson, B., Gustafsson, B., 1985, MNRAS, 212, 497. 

Berdyugina, S. V., Savanov, I. S, 1992, SvA, 36, 425. 

Bonnell, J. T., Bell, R. A., 1993, MNRAS, 264, 334. 

Burrows, A., Ram, S.R., Bernath, P. et al. 2002, ApJ., 577, 986. 

Burrows, A., Volobuyev, M. 2003, ApJ, 583, 985. 

Chabrier, G., Baraffe, I., Allard, F., and Hauschildt, P., 2000, ApJ, 542, L122. 
Craybeal J.D. 1988. Molecular Specroscopy, McGraw-Hill Book Company. 
Engvold, O., Wohl, H, Brault, J.W., 1980, A&A Suppl. Ser., 42, 209. 
Gaydon, A.G., Pearse, R.V.B. 1937. Nature, 140, 110. 



Granovsky, A. A. 2004, http://quantum-2.chem.msu.ru/gran/gamess/index.html 



20 Ya. V. Pavlenko et al. 

Gurvitz, L. V., Weitz, I. V., Medvedev, V. A. 1989, Thermodynamic properties of individual 

substances. Moscow. Science 

Halfen, D.T., Zirius, L.M. 2004. 611, L65. 



Hill, C. 2007 http://www.physchem.ox.ac.uk/liill/researcli/MgH/index.html 
Jones, H.R.A., Tsuji, T. 1997, ApJ, 480, 39. 

Kirkpatrick, J. D., Allard, F., Bida, T., Zuckerman, B., Becklin, E. E., Chabrier, G.; Baraffe, 
I. 1999, ApJ, 519, 834, 

Kleman, B., Uhler, U., 1959, Can. J. Phys. 37, 537. 

Kuz'menko, N., Kuznetzova, L.A., Kuz'yakov, Y.Y, 1984, Franc-Condon factors for two- 
atomic molecules, Moscow University Publ., 1. 

Kuznezova, L.A., Kuzmenko. N.E., Kuz'yakov, Y.Y., PLastinin, Y.A. 1989, Two atomic 

molecules transition probabilities, Moscow, Nauka, 315 p. 

Leggett, S. K., Allard, F., Geballe, T.R., at al., 2001, ApJ, 548, 908. 

Lemoine, B., Demuynck, C.,Destombes, J.L., Davies, P. B., 1988, J. Chem. Phys, 89, 673 

Lipus, K., Bachem, E, Urban, W. 1991, Mol. Phys., 73,1041. 

Magazzu, A., Martin, E. L., Rebolo, R., 1993, ApJ, 404, L17 

Martin , E. L., Delfosse, X., Basri, G., Goldman, B., Forveille, T., Zapatero Osorio, M.R. 
1999, AJ, 118, 2466. 

Naersisian, S. E., Shavrina, A. V., Yaremchuk, A. A., 1989, Astrofisika, 30, 147. 

NichoUs, R.W. 1981, ApJSS, 47, 299. 

Pavlenko Ya.V. 1997, Odessa Astron. Publ., 10, 76. 

Pavlenko, Ya. V. 1997a, Astron. Astrophys. Sci., 253, 43. 

Pavlenko Ya.V. 1998, Atronom. Report, 42, 787. 

Pavlenko Ya.V. 1998a, Atronom. Report, 42, 787. 

Pavlenko, Y. V. 1999, Astron. Rep., 43, 748. 

Pavlenko, Ya.V. 2000, Astron. Rept., 44, 219 

Pavlenko, Ya.V. 2002, Astron. Rept., 46, 567 

Pavlenko, Ya.V., Zhukovska, S., Volobuyev, M. 2007, Astron. Rept., 51, 282.. 

Ram, R.S., Bernath, P,F., 1995, J. Molec. Specrosc, 172, 91. 

Ram, R.S., Jarman, C.N., Bernath, P,F., 1993, J. Molec. Specrosc, 161, 445. 

Rebolo, R., Martin, E. L., Basri, G., Marcy, G. W., Zapatero-Osorio, M. R. 1996, ApJ, 

469, 53. 

Saumon, D., Hubbard, W.B., Burrows, A. et al. 1996, ApJ, 460, 993 



Bands of CrD and MgD and the "deuterium test" 21 

Schadee, A. 1968, ApJ, 151, 239. 

Shaycstch, A., Henderson, R. D. E., Le Roy, R. J., Bernath, P. F. (in preparation). 
Shin, S., Bruch, D.J., Morse, M.D. 2005, ApJ, 619, 407. 
Sinha, K., Shukla, D. S., Tripathi, M. B. 1979, BASI, 7, 38. 
Sinha, K., Joshi, G. C, 1982, BASI, 10, 329. 

Skory, S., Week, R F., Stancil, P. C, Kirby, K. 2003, ApJS, 148, 599. 

2003, ApJ, 584, 459 

Tsuji, T. 1973, A&A, 23, 411 

Tsuji, T., Ohnaka, K., Aoki, W., Nakajima, T. 1996, A&A, 308, L29. 
Tsymbal, V.V. 1977, Tables of Franck-Condon factors with account the rotational- 
vibrational interactions for astrophysically important molecules. I. Titanium Oxide, Dep. 
No. 246-77. 

Zu-Geng Wang, Hui-Rong Xia 1996, Molecular and Laser Spectroscopy, Springer Verlag. 

Week, P. F., Schweitzer, A., Stancil, P. C., et al. 2003, ApJ, 584, 459. 

Week, P. F., Stancil, P. G., Kirby, K, 2003a, ApJ, 582, 1263. 

Week, P. F., Schweitzer, A., Stancil, P. G., et al. 2003b, ApJ, 582, 1059. 

Wyller, A. A., 1961, ApJ, 134, 805 

Ziurys, L. M., Barclay, W. L. Jr., Anderson, M. A. 1993, ApJ, 402L, 21. 
Zhu, G., Babb, J.F., and Dalgarno, A. 2006, Phys. Rev. A, 73, 012506. 



, K(D+MgH / H+MgD) 

K(D+CrH / H+CrD) 

0.2 L 



500 1000 1500 2000 2500 3000 3500 4000 4500 5000 

T[K] 



